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Damped harmonic oscillator interpretation of the soft-state power 
spectra of Cyg X-l 
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ABSTRACT 

We develop a model of an accretion disc in which the variability induced at a given radius 
is governed by a damped harmonic oscillator at the corresponding epicyclic frequency. That 
variability induces both linear and non-linear responses in the locally emitted radiation. The 
total observed variability of a source is the sum of these contributions over the disc radius 
weighted by the energy dissipation rate at each radius. It is shown that this simple model, 
which effectively has only three parameters including the normalization, can explain the range 
of the power spectra observed from Cyg X-l in the soft state. Although a degeneracy between 
the black hole mass and the strength of the damping does not allow a unique determination of 
the mass, we can still constrain it to < (16-2O)M0. We also show that our model preserves 
the observed linear rms-flux relationship even in the presence of the non-linear flux response. 
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1 INTRODUCTION 

Black-hole binaries ar e primarily observed to be i n two main spec- 
tral states (see, e.g., [zdziarski & Gierlinski 2004| for a review on 
the spectral states and the related variability). In the hard/low state, 
the spectrum can be approximately described as a hard power-law, 
with the typical photon index of T » 1.7, and a high energy cut- 
off above ~ 100 keV (Gierli nski etafl 1 1997b . and only a weak 
soft X-ray component. The likely physical origin of that spectrum 
is thermal Comptonization of some soft photons by a hot plasma 
with kT ~ 10 2 keV (e.g.. iGierlifiski et all \l99% I Frontera et al.1 
l200ll) . The hot plasma may be identified as an inner hot accre- 
tion disc surrounded by a truncated outer cold disc, which provides 
seed photons for Comptonization dShapiro. Lightman & Eardlevl 



1 19761 : lEsin. McClintock & Naravanl 1 1991 IZdziarski et alj|2002h 
On the other hand, in the soft/high state, the spectrum consists of 
a soft blackbody-like component, which dominates the luminosity, 
and a hard X-ray power-law tail with a photon index T ss 2.5- 
3 with no detectable high energy cutoff. While the soft compo- 
nent is consistent with the emission of a cool accretion disc ex- 
tending close to the innermost stable orbit, the nature of the hard 
power-law tail is less clear. Detailed modeling of spectra and 
long-term variability shows that it is well represented by Comp- 
tonization of soft photons by a hybrid plasma, i.e., containing 
both thermal and non - thermal electrons Jpoutanen & Coppi|[T998l : 
IGierlifiski etafl Il999l : IZdziarski et al.l l200lL l2002h 7 This plasma 
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may exist in the form of active, time-varying, regions above the 
cold disc faaardt & Maraschilll993l : |Poutanen & Fabianlll999h . In 
addition, during transitions between those two states, black-hole bi- 
naries also show the state called intermediate or very high. In each 
state, there are significant spectral and flux variations on a variety 
of time-scales. 

The temporal variability of a black hole system is different in 
the hard and soft spectral states. During the hard state, the power 
spectrum, P, is approximately a broken power law with the breaks 
at the frequencies of / ~ 0. 1 Hz and / ~ 5 Hz dNowak et all 19991: 
Gilfa nov. Churazov & Revnivtsev|[l999h . whereas the power spec- 
trum in the soft state is approxim ately a power law with a high fre- 
quenc y cutoff above / ~ 10 Hz dChurazov. Gilfanov & Revnivtsevl 
l200ll) . In addition, discrete features, referred to as quasi-periodic 
oscillations (QPOs), are often observed. Attempts have been made 
to identify the centroid frequencies of these QPOs with frequencies 
at spe cific radii of the disc, w hich could either be the truncation ra- 
dius (IGiannios & Spruitl2004l) or the radius of a resonance between 
the r adial and vertical epicyclic frequencies fAbramowi cz et all 
2003). These intriguing interpretations do not generally address the 
origin of the continuum variability. Detailed but empirical fitting of 
the power spectra using Lorentzian functions reveal that the con- 
tinuum variability is complex in nature. Often three or more com- 
ponents (with a total of nine or more param eters) are required to 
adequately describe the spectra dNowakll2000li . However, it is pos- 
sible that this apparent complexity appears because the Lorentzian 
function may not be an adequate description of the variability. In 
the framework of the active flare model, the continuum variabil- 
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ity can be explained by the rise and f all time scales of the flare s 
and the frequency at which they occur dPoutanen & Fa bian 1999). 
However, this interpretation does not make a direct connection with 
the characteristic dynamic frequencies of the accretion disc. 

The accretion process in a disc occurs via angular momentum 
loss through turbulence, which nature and origin remain uncertain . 
In the standard accretion disc model (Sh akura & Sunvaevlll973h . 
this uncertainty is parametrized by approximating the turbulent be- 
havior as an effective viscosity providing a stress proportional to 
the local pressure. While such a prescription may indeed provide 
a reasonable approximation to the time-averaged disc, the effect of 
turbulenc e on time-dependen t properties of the disc is much more 
uncertain. lLvubarskii|]T997h showed that if the viscosity parame- 
ter varies at each disc radius with a local characteristic frequency, 
the accretion rate in the inner disc fluctuates with the flicker-noise 
behaviour, i.e., with P oc 1//. This provides a natural explanation 
for the observed power-law of the soft-state of black hole systems. 
Also, it explains the X-ray variability observed mostly at low fre- 
quencies in spite of the characteristic frequencies of the inner disc, 
where the X-rays originate, being much higher. Such a generic 
model where variability propagates from larger radii to smaller 
ones can also explain t he observed frequency-dependent photon- 
energy lags in Cyg X-l dMisrakOOrilkotov. Churazov & Gilfanovl 
l200ll) . It is also consistent with th e observed linear rms-flux rela- 
tionship dUtflev & McHardy||200ll) . 

In these models, a high frequency break in the power 
spectrum is likely to correspond to a characteristic frequency 
at the innermost radius. Since the break obser ved in the soft 
state is typically at ~ 10 Hz, ILvubarskiil d 19971) suggested that 
this corresponds to the viscous t i me sc ale of the inner disc. 
iTitarchukTs haposhni kov & Arefievl d2007l) have generalized this 
model by invoking a power-law variation of the viscosity parame- 
ter with radius, and used a specific form for the local perturbations. 
They show that the observed power spectra of black hole systems 
can be explained by a model with two components, corresponding 

to the outer disc and an inner hot region. 

A crucial element of such models, as noted by ILvubarskiil 

(1997), is that the time scale for the variations should be as long 
as the viscous one in order to explain the observed high frequency 
break. On the other hand, the main charact eristic time scale is the 
epicyclic one, which is significantly shorter. King et al] d2004l) note 
that even the characteristic magnetic dynamo time scale is also 
expected to be significantly smaller than the viscous one. If the 
dissipation is driven by turbulence induced by magnetic instabili- 
ties, it is that dyn amo time scale w hich should dominate t he ob- 
served variability. iKing et alj J2004I) and iMaver & Pringld J2006I) 
provide a solution to this problem by suggesting that large scale 
outflows can only occur when weak random events in different disc 
annuli by chance give rise to a coherent poloidal field. The time 
scale for that is much longer than the local dynamo one. It is then 
these large scale outflows that govern the observed slow variability. 
Those authors have developed a specific disc model incorporating 
the stochastic dynamo action as well as jet/wind outflows. Their nu- 
merical results suggest that it may also explain observed variability 
of black hole systems. 

In all these models, the variability is associated with the un- 
known turbulence and not with the characteristic dynamic time 
scales of the disc. On the other hand, the simplest variability model 
would associate it with epicyclic oscillation at each radius, with 
the net variability due to a superposition of these oscillations. 
In this approach, even if variability at large r adii induces som e 
low frequency fluctuations in the inner regions dLvubarskii|[l997h . 



the dominant local power should correspond to the epicyclic time 
scales. However, such a natural model in its simple version would 
predict a much higher break frequency than that observed. 

In this work, we propose that these natural oscillations do in- 
deed occur, but a damping, viscous-like, effect, suppresses the high 
frequencies. A likely reason for the viscous damping is the same 
turbulence that causes the viscous dissipation. Thus, contrary to 
the models mentioned above where the observed variability is pri- 
marily due to variations in the turbulent viscosity and/or dynamo 
effects, the variability in our model is due to the natural dynamic 
oscillations of a disc but damped by viscous-like effects (perhaps 
also the turbulence). Based on this idea, we have developed a sim- 
ple framework that explains the power spectra of black-hole bi- 
naries in the soft state. We limit the analysis to this state, where, 
as mentioned earlier, an optically-thick disc extends nearly to the 
last stable orbit. Although most of the observed variability in the 
soft state occurs in the high energy photons, and hence it should 
be associated with active regions or flares, it is expected that the 
underlying disc is still the main driver of the variability (e.g., via 
magnetic buoyancy). 

We assume that the disc variability at each radius is governed 
by a damped harmonic oscillator at the corresponding epicyclic fre- 
quency. This dynamic variation induces both linear and non-linear 
responses in the emitted radiation. The total observed variability of 
a source is then a sum of these contributions, weighted by the rela- 
tive contribution of the radiation from a given radius. We find that 
the predicted power spectrum is completely determined by only 
three effective parameters including normalization, and it can ex- 
plain the different spectral shapes observed for Cyg X-l in the soft 
state. Whereas a degeneracy between the black hole mass and the 
damping factor does not allow the determination of the mass of the 
system, the model constrains it to < (16-20)Mq, which is consis- 
tent with other mass measurements. 

In Section [2] we describe the model. In Section [3] we com- 
pare the predicted power spectra with observations from Rossi X- 
ray Timing Explorer (RXTE). In Sectionf4] we discuss implications 
of our results for the relationship between the rms and the flux. In 
Section|5] we briefly discuss and summarize the results. 



2 THE MODEL 

The variability of a dynamic quantity, X(t), at each radii of the ac- 
cretion disc is assumed to be governed by a damped harmonic equa- 
tion, 



m + fi(r)x(t)+Mr)m = C(t), 



(1) 



where x(t) = [X{t) - X ]/X is the normalized fractional varia- 
tion, with .Xo being an average of X(t). Then C(i) is a stochastic 
driving term and /e(^) is the epi cyclic frequency, wh ich for the 
Schwarzschild metric equals (e.g. JStella & Vietf 



c y, wm c 
3l999h - 



1/2 



(2) 



where r g = GM/c 2 and M is the mass of the black hole. Then 
fa is an unknown damping frequency, which we parameterize as 
fo( r ) - Af E (r) (A = corresponds to no damping). While we do 
not specify the nature of the dynamic normalized variable, x(t), we 
assume that it induces a linear response departing from the average 
photon flux, So, such that s(t) = [S(t) - So] /So oc x(t). Thus, the 
observed linear compone nt of the power spectrum corresponding 
to an annular radius, r, is dLandau & Lifshitzll 19691) . 
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Pl(f) = Mf)\ 2 Wftf oc 



|C(/)I 2 



(/ E 2 -/ 2 ) 2 +^l 2 / E 2 / 2 ' 



(3) 



where .?(/), x(/) and C(f) are the Fourier transforms of s(f), x(t) 
and C(f), respectively. It is assumed that in the frequency range of 
interest, the driving fluctuations, C(t), are dominated by a 1// noise 
such that |C(/)| 2 oc f~' valid for / greater than a low frequency 
cutoff, yj„i n . Thus, 

a- / 



where the normalization, 



r'd/ 



(4) 



(5) 



lf mm (fi-f 2 ) 2 + * 2 flf 2 

makes the fractional rms of the observed variability equal to a. If 
/min ^ / (which we assume hereafter), / becomes, 



Hh,X) * / E 4 pn(A//«dn) + B(-i)] . 



where 

f [b(A)/4] [n + 2 arctan 6(A)] , 1 < 2; 
g( /i )= -1/2, /f = 2; 

and 

2 -A 2 



A>2. 



b(A). 



a^\a-a 2 \ 



(6) 



(7) 



(8) 



Thus, the normalization of p\(f) depends only logarithmically on 
/min> while the spectral shape is independent of it. 

The total observed power spectrum, P\(f), is then a weighted 
sum of the contributions from different radii, 



Hf) : 



I™Pi(/)G(r)27trdr 
£° G(r)2nrdr 



(9) 



where G(r) is a weight function describing the relative contribu- 
tion of each radius to the total variability, p\ is also dependent on 
r via the epicyclic frequency, equation (O, and on A. We assume 
hereafter r m = 6r g . In writing equation {5), it has been implicitly 
assumed that the driving fluctuations at different radii, C(t), are in- 
coherent. This assumptions requires that r coh <k r, where r coh is 
the radial coherence length of the underlying variability, which is 
justified if, e.g., r CO h is of the order of the scale height of the disc. 
Most of the variability observed in black-hole sys tems in the soft 
state is associated with the high energy component dChurazov et al.l 
1200 lh . which probably arises from hot active regions on top of the 
cold disc. The radial extent of these active regions is not clear, and, 
in particular, the energy dissipated in them as a function of radius 
remains unknown. Thus, it is difficult to obtain the form of the 
weight function, G(r), from first principles. Here we mak e a simple 
assum ption (also used for modeling disc rms spectra inlZdziarski 
l2005t) that G(r) is proportional to the local energy dissipated in the 
cold disc, i.e., 



G(r) oc r 



6r„ 



1/2 



(10) 



The observed spectrum, P\(f), depends on the normalization 
factor, a, the damping factor, A, and, through / E , on the black-hole 
mass, M. Of these, a is the normalization, while variations in M 
shift the spectrum in frequency. Hence, A is the single parame- 
ter that describes the spectral shape in this linear model. In Fig. 
[]2a), we show Pi for several different values of A and M = IOMq, 




0.1 1 10 100 

Frequency (Hz) 

Figure 1. The power per logarithm of frequency, fP(f), for M = IOMq. 
The curves from right to left are for A = 0.1, 0.3, 1, 3, 10, 30, 100, for (a) 
the linear power, P\, and (b), the non-linear one, P n \. The normalization of 
each curve corresponds to unity at (a) /'Pi(0.01 Hz), (b) the peak of fP n \. 



,/min = 10 4 Hz. At small values of A <k 1, the oscillator is almost 
undamped, and the spectrum has a narrow spike at 71(10Mq/M) 
Hz. This corresponds to the maximum epicyclic frequency, fc m , 
equation {2), which is achieved at r - 8r g , very close to the max- 
imum of the gravitational energy dissipation rate, G(r), equation 
d 1 0b . In the time domain, we can define the corresponding time 
scale as (27t/ Em ) _1 ^ 2(M/10Mq) ms. Observationally, it roughly 
corresponds to the shortest e-folding time scales seen yet (which 
happened in the soft state) from C yg X-l dGierlinski & Zdziarskil 
2003: lZdziarski & Gierlifiskill2004l) . 

For / <K / Em , P] oc l/f, which is a consequence of the form 
of driving variability of C(f) oc l/f. As A increases, the system 
becomes more damped, and the spectrum shows a cutoff at a fre- 
quency <K / E m . For large values of A > 1, an increase in A is nearly 
equivalent to a scale shift of P\(f) in frequency. Such a frequency 
shift also occurs when the mass of the black hole is increased. This 
results in a degeneracy between the parameters, M and A, when the 
model is fitted to observational data (see Section[3)- Therefore, we 
can effectively determine only A/M. 

The response of the observed flux to the variability, x(t), may 
also have a significant non-linear component, which in the leading 
order would be s n ](t) oc x 2 (t) - (x 2 (t)). The Fourier transform of j„i 
can be expressed as, 



j(/) ' 



£*nxf-rw. 



(id 



Expressing x(f) = -y^e*-" and changing the variable to A/ = 
f/2 — /', the non-linear term can be rewritten as, 

5„,(/) oc J Vpi(//2 - A/)p,(//2 + A/)e* //2+A/) - w/2 - A/) dA/.(12) 

Since the driving term C(f) is stochastic, the phase, (f>(f), will be 
largely incoherent. Thus, the above integral will average to zero 
except when 2A/ < A/ coh , where A/ coll is the frequency range over 
which the phase of C(f) remains coherent. For small A/ co h <k /, 
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S„i K Pi(//2)A/ coh , and the non-linear power spectrum can be ex- 
pressed as 



P„l(/) ' 



We,-*) 



fpt(f/2), Wb,2): 



I 



fp;(f/2)df, (13) 



where A/ coh oc / has been assumed, I nl is the normalization, and /? 
is the fractional rms. 

The observed non-linear component of the power spectrum, 
P„i, is then the weighted average of p n \ over the energy dissipa- 
tion distribution, G(r), analogously to equation l[9}. Apart for the 
normalization factor, f3, P nl depends on A and M. Similar to the lin- 
ear case, variation in M results in an overall frequency shift, and 
the shape of the spectrum is determined only by A. In Fig. |TJb), 
the computed spectra, P n \, are shown for different values of A and 
M = IOMq. Similar to the linear component, variation of A when 
A > 1 results only in a frequency shift, which is degenerate with 
variation of M. 

The total observed spectrum will, in general, be a combination 
of the linear, Pi, and non-linear, P nl , components, with the total 
fractional rms of -\/(a 2 +fP). Apart from this normalization factor, 
the spectrum depends on A, the black hole mass, M, and the power 
ratio between the non-linear and linear components, q nl = /3/a. 



3 FITS TO DATA 

The Proportional Counter Array (PCA) on board RXTE has repeat- 
edly observed the Cyg X-l. During some of the observations the 
system was in the soft state, and a criterion to select soft state 
data can be obtained by defining a soft colour, C s , as the ratio 
of th e energy fluxes in the en ergy bands 4.0-6.4 keV to 3.0-4.0 
keV jPone & Gierlihski 2003). The system exhibited a range of 
0.7 < C s < 2.4, and in this work we define the system to be in 
the soft state when C s < 1.2. We have selected 1 1 pointed observa- 
tions (during 1996-2004) satisfying this criterion and covering the 
range of 0.7 < C 8 < 1.2 approximately uniformly (M. Gierlihski, 
personal communication). TableQ]gives the IDs of the selected ob- 
servations. Rebinned power spectra were calculated from the light 
curves using standard techniques. A 5 per cent systematic error was 
added in quadrature to the statistical error of the spectra. 

The spectra were then fitted to the damped harmonic oscillator 
model for / > 0.1 Hz. The resulting best fit parameters and their 
90 per cent uncertainties (i.e., for A^ 2 = +2.71) are given in Table 
Q] Four representative spectra and their fits are shown in Fig. [2] 
As mentioned in Section[2] there is a near degeneracy between the 
black hole mass and the damping factor, A, in all cases. Therefore, 
the fits in Fig. [2] are for M = IOMq. Generally, the best fits are 
obtained for low masses, and we can obtain only upper limits, M u , 
using the criterion of Ax 2 = +2.71 with respect the overall best fit 
for each data set. As we show in TableQ] we find M u in the range of 
(16-36)Mq, though M„ < 20Mq in most cases. We illustrate this 
degeneracy in Fig. [3] where we plot the variations of the reduced 
X 1 and of A vs. M for one of the two data sets yielding the lowest 
M u ss 16Mq. Thus, the limit of M < 16Mq is satisfied by all of the 
data, and we can consider it as our resulting overall limit. On the 
other hand, M < 20Mq is satisfied by most of our data. 

The mass of the black hole in Cyg X-l determ ined using var- 
ious o ther techniques remains relatively uncertain. iGies & Boltonl 
( 1986) foun d M ^ 1 6 ± 5M ^, consistent with the later determi- 
natio n by Izioikowskil d2005l) of 13 < M/Mq < 29. On the other 
hand, iHerrero et alj 1 19951) . based on modeling of the atmosphere 




io- a io-« 10° 10' to 2 io- 2 to- 1 10° 10' 10 2 
Frequency (Hz) 

Figure 2. The observed power spectra and best fit models (solid curves; 
fitted at / > 0. 1 Hz) for four representative data sets. The dotted and dashed 
curves give the linear and non-linear components, respectively. 



Table 1. Results of the fits to soft-state power spectra of Cyg X-l. The 
values and errors of A, are given for M = IOMq; those for fractional rms 
and q n i are approximately independent of M. 



Observation ID 



rms (%) 



<?nl 



Xv 



MJMq 



10512- 
10512- 
50110- 
50110- 
50110- 
60089- 
60090- 
60090- 
60090- 
60090- 
70015- 



01-09-00 
01-09-01 
01-52-00 
01-44-01 
01-53-00 

03- 01-01 
01-07-01 
01-11-00 
01-12-00 
01-15-01 

04- 01-00 



5 22 

5 02+^ 
-Q-2Q 

5 [7+0.89 

<: 1 0+8:^8 
4 4<j+0.82 

6 31+o.W 

D - J1 -0.49 
=; 41+0.69 

497^ 

^■ y/ -0.65 
4 63+ - 24 
c Q4+0.43 
J - -0.40 



24!] 
32+| 
36+1 
27+! 
28 + ] 
35+! 



24+ 
36_ 
33 + j 
27+| 
23+! 



+1 



46+ 004 
0.4 1+ 003 

13+^ 
-Q-P? 

57+™ 
-Q-Q3 

\5 +m 
u - -0.03 

22+^+ 
U - 4S -0.04 

13+ - 03 

15+™ 
u -0.03 
n 44+0.03 
U - 4 ^-0.03 

61+ 06 

-0.05 



1.27 
1.10 
0.86 
0.89 
1.13 
0.82 
1.30 
1.16 
0.84 
0.77 
1.27 



20 
18 
16 
18 
18 
24 
28 
16 
18 
22 
36 



of the companion star, obtained M ^ 10 ± 5Mq. Still, those deter- 
minations are consistent with our limit of M < (16-20)Mq. 

As seen in Fig. [3] the fit in the shown case yields the damping 
factor of A > 3. On the other hand, we can use the independent 
limits on M to cons train A from above. In any c ase, M > 5Mp 
( Herre ro et al. I GUI), or, more likely, > 6.5Mq ( Gierlihski et al. 
1999). The latter implies A < 8. These constraints are consistent 
with TableQ] which shows A ~ 5 ± 1 for M = IOMq, taking into 
account that our fits constrain mostly A/M (Section|2j. 



4 THE RMS-FLUX RELATIONSHIP 

Above, we have shown that our model fits well the soft-state power 
spectra. Another major observable describing variability of accret- 
ing sources is the flux rms as a function of the flux, which has been 
shown to be linear in accreting black-hole and neutron-star systems 
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5 DISCUSSION AND CONCLUSIONS 



5 i 1 , 



* * i , 



M (M s ) 

Figure 3. The variations of x 2 and the best fit damping factor, A, vs. the 
black-hole mass for the observation 501 10-01-52-00 (see Fig. [2), which is 
one of the two data sets yielding the lowest M a . At M > M u = 16Mq, 
Ax 2 > 2.71 with respect to the overall^ 2 minimum, as shown by the solid 
horizontal line in the bottom panel. Hence, I6M0 is an upper limit for M. 
The corresponding constraint for the damping factor is A > 3. 



dUttlev & McHardvl200ll) and AGNs (e.g. jMarkowitz et alJ2007h . 
In particular, the fr actional rms ^ const ant in both the hard and soft 
states of Cyg X-l ddeissner et alj|2004h . 

We point out here that, in our formalism, this linearity is an 
intrinsic property of the dynamic quantity governing the variability, 
x(f), but it is not related to the relationship between the flux and x(t). 
The flux variance over a time interval, T, is 



- [S(t)-s Ydt, 

1 Jo 



(14) 



or the corresponding sum over some time steps, Ar. If the normal- 
ized flux is s(t) = a [x"(t) - (x"(t))], the flux variance becomes, 



£i f 7 

T Jo 



[x"(t)-(x"(t))fdt, 



(15) 



where n = 1 in the linear case, and n > 2 in the non-linear one. 
Then, if the variance of x(t) [given by the integral in equation U5i 
for n = 1] is approximately the same in different time intervals, 
T, i.e., cr(X) oc X in each time interval, it is easy to see that 
<t(S) oc So at any n, with only the proportionality constant chang- 
ing. As an illustrative example, let's take a single harmonic of x(f), 

1. e., x(i) = sin/?, for which a 2 = and Sq/8 for n = 1 and 

2, respectively, i.e., the rms-flux relationship remains linear even 
when the flux response to the underlying variability is non-linear. 

In our model, the power is supplied to the disc via the driving 
term, C(t). As x(t) is linearly related to C(f), the observed rms-flux 
linearity can occur if the standard deviation of C(t) is proportional 
to its local average. Since we do not specify the nature of C{t) [ex- 
cept for assuming C(f) oc 1//], our model does not imply a partic- 
ular rms-flux relation, but it is consistent with it being linear. 



We have developed a relatively simple accretion disc variability 
model, in which we assumed that the local dynamic variability at 
each radius is governed by a damped harmonic oscillator equation 
at the corresponding epicyclic frequency. This variability induces 
both linear and non-linear (quadratic in our model) responses in 
the emitted photons. The observed variability is then the sum of 
those from different radii weighted by the rate of local energy dis- 
sipation. The power spectrum predicted by our model is charac- 
terized by four parameters, which are the total fractional rms, the 
ratio of the non-linear power to the linear one, q ut , the mass of the 
black hole, M, and the damping factor, A. We find that M and A are 
closely correlated, and we can determine only A/M. Thus, the pre- 
dicted spectrum is a function of only three effective parameters. We 
find we can well explain the different types of the power spectra ob- 
served from Cyg X- 1 in the soft state with this model. Although the 
degeneracy in A and M does not allow an unambiguous determina- 
tion of M, we are able to obtain an upper limit of M < (16-20)Mq. 
Remarkably, this limit is close to the range of other mass measure- 
ments of Cyg X-l. The corresponding constraint on the damping 
factor is 3 < A < 8. 

This model should, in the future, be made consistent with a ra- 
diative model for the variability. Whereas we here assume that the 
dynamic disc variability induces linear and non-linear responses in 
the emitted radiation, a radiative model should identify the parame- 
ters of the emitting medium (e.g., the optical depth, the plasma tem- 
perature, etc.) that drive the variability. It should predict the rms vs. 
the photon energy and time lags between different energies. One 
important issue we need to understand within such framework is 
the nature of the non-linear response of the plasma, which inher- 
ent presence in the soft state of Cyg X-l we have found from the 
observations. On the theoretical side, we can state that its relative 
strength, q n \, depends on the frequency range, A/ CO h, over which 
the stochastic driver C(t) is coherent. If the nature of the stochas- 
tic driver does not change, variations in the strength of the non- 
linear component imply that the non-linear response of the plasma 
varies. This predicts q nl to be correlated with the parameters of the 
hot plasma determining its emission. However, quantitative predic- 
tions can be made only when a comprehensive radiative model for 
the variability is available. 

An important constraint on C(t) is that it should preserve the 
observed linear relationship between the rms and the flux. As we 
have shown, this implies that the variance of C(t) remains approx- 
imately constant over time. On the other hand, we have shown that 
even a non-linear response of the flux to the underlying dynamic 
variability preserves the linearity of the rms-flux relation. 
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